Groundstates of 5'[/(2)-Symmetric Confined Bose Gas: Trap for a Schrodinger Cat 
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Conservation of the total isotopic spin S of a two-component Bose gas — like 87 Rb — has a 
dramatic impact on the structure of the ground state. In the case when S is much smaller than 
the total number of particles N, the condensation of each of the two components occurs into at 
least two single-particle modes. The quantum wavefunction of such a groundstate is a Schrodinger 
Cat — a superposition of the phase separated classical condensates, the most "probable" state in 
the superposition corresponding to the classical groundstate in the sector of given N and S. 

PACS numbers: 03.75.Fi 



Symmetry has profound implications on the ground- 
state structure of the degenerate superfiuids and super- 
conducting systems. A remarkable example to that is 
superfluid 3 He, which exhibits a variety of textures |Q. 
These ground states are essentially classical — quantum 
fluctuations play no role at large length scales. Accord- 
ingly, the most important low energy physics can be well 
described within the mean-field (MF) approach relying 
on broken gauge symmetry. 

Recent successes in trapping and manipulating atomic 
Bose-Einstein condensate (BEC) have initiated a 

great interest to the physics beyond the MF in these sys- 
tems. A well-known example of the non-MF behavior is 
the phase diffusion effect 

A remarkable case of non-MF ground state was found 
in the spin s = 1 sodium vapor |||(| . This state is charac- 
terized by anomalous quantum fluctuations of the com- 
ponents with different projections of the total spin. 

In this Letter, we address the case of a two-component 
SJ7(2)-symmetric Bose gas — like 87 Rb j7|, that is, the 
case of (pseudo-) spin- 1/2 bosons. We show that, while 
sharing some common features with the s = 1 bosons 
H||, the s = 1/2 system exhibits unique properties ||. 
We pose a question of the formation of the ground state 
in a naturally arising situation of S <C N. We find that 
the resulting ground state is built on at least two macro- 
scopically populated lowest one-particle modes, and is 
a quantum superposition of the phase separated BECs, 
characterized by large fluctuations of the densities of the 
components, the total density being fixed. 

We start from the S[/(2)-symirietric Hamiltonian 



(i) 



where the interaction constant g = Aita/m is repre- 
sented (in units K = 1 ) by the scattering length a and 
mass m; the single particle Hamiltonian Hi con- 
sists of the kinetic energy and the trapping potential, 
H 1 = -(l/2m)V 2 + f7(x); the Bose-field operators 1^ 
have two components a =f, J, in pseudo-spin notations, 
and the summation is performed over repeated pseudo- 
spin (Greek) indices. Due to its symmetry, this Hamilto- 



nian conserves the total isospin S = (S x , S y , S z ) given 
as ( S± = S x ± iS y , S+ = S_) 



±JdK(9\9i-*l9 l ),S-=JdK*l* r . (2) 



S '=2 



These operators obey the standard SU(2) algebra 
[S Z ,S+] = S +1 [S„S-] = -S-, [S+,S^} = 2S Z . 

The absolute minimum of energy is achieved when all 
atoms condense to a lowest one-particle (Hartree) state 
ipo. Such a state corresponds to the maximum possible 
value of the isospin S max — N/2 = (JV-j- + -/VjJ/2, where 

, N± stand for the total numbers of atoms of each com- 
ponent. Thus, it follows that a very nontrivial situation 
will occur when, while conserving S, the condensation 
proceeds from the thermal cloud which contains almost 
the same amount of different components N a with no 
inter- component correlations. Indeed, the ensemble mean 
of the total isospin (|J) in Boltzmann gas is (S) = 0, 
and the mean of S 2 = S 2 Z + (S+S- + SS+) /2 yields 
(S 2 ) - N, as long as (*{(x)* i (x')) = and the di- 
agonal correlators decay on thermal length. Similarly, 
(S 2 ) ~ N, when two spatially separated degenerate com- 
ponents are mixed together (see below). Thus, the re- 
sulting ensemble of the ground states is characterized by 
typical values < S ~ y/N. 

Since condensation of both components into one single- 
particle state <po is incompatible with S <C N, the result- 
ing many-body state is characterized by macroscopic oc- 
cupation of at least two lowest one-particle states. This 
contrasts the cases s = 1/2 and s = 1 

We analyze the simplest case of the strong quasi- ID 
trap, where the energy difference e = e\ — Eq between 
the lowest single particle levels tpifi, respectively, in the 
trapping potential U is much larger than a typical in- 
teraction energy per particle \i (chemical potential at 
T = 0). The inequality /j/e<1 allows one to essen- 
tially simplify the description by introducing a projected 
Hamiltonian acting in a truncated Hilbert space contain- 
ing only two single-particle modes, ipo and ip\: 



*cr = a brT ip b 
6=0,1 



(3) 
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where a ba annihilates a particle at the level b in the 
spin state a. With the same accuracy, one may ne- 
glect in the projected Hamiltonian all the terms that do 
not conserve the total occupations of each of the two 
single-particle levels. These terms involve a large energy 
difference ~ e ^> n, and, thus, result only in higher- 
order — in parameter \i/e — perturbative renormalization 
of the parameters of the effective (two-modes & two- 
colors) Hamiltonian 



o=0,l 

+ i Y, IabN a (N b ~S a , b ) + 2I Q1 S S u (4) 



,6=0,1 



where I ab = g J dx\ip a \ 2 \ip b \ 2 and N b = J2a N bcn S bz = 
(A hT - N H )/2, S bx = (a^ahi + a\ i a b {)/2, S by = 
(aLatj. — ahab-\)/2i; the operators N aa = ol a a ba and 
S a = (S ax , S ay , S az ) represent the numbers of atoms and 
the components of the pseudo-spin on the ath level, re- 
spectively. The quantities N a and S 2 are related to each 
other by the identity Si = N a (l + N a /2)/2, implying 
that S a = N a /2. 

A standard choice of eigen numbers for the Hamilto- 
nian (Q) comes from the theory of two spins: the total 
spin S, its z-projection, M, and the spins So and Si; with 
the constraints |S - S x \ < S < S + Su -S < M < S. 
In this nomenclature, the spectrum of (|]) is 

EsmSqSi — 2eo('5o + Si) + 2sSi + ^ I ab S a (2S b — S a> b) 

a . 6 

+/oi[S(S + 1) - (S + 1)S - (Si + l)5i] . (5) 

The distribution over the eigennumbers So, Si, S, M, 
which will be formed in the end of the cooling pro- 
cesses, depends essentially on the following: While S, 
M and N are the eigennumbers of H, eq. ([!]), the spins 
So, Si are not the eigennumbers of the original Hamil- 
tonian (|l]). Thus, the formation of particular low energy 
states, amenable to treatment by the reduced Hamilto- 
nian (§) where S ,i are now the quantum numbers, 
will be controlled by the microcanonical distribution with 
fixed values of S, M, N. For our purposes, we do not 
need to adopt any particular distribution for S, M and 
N. We assume that their values are obtained (after the 
cooling is over) by a non-destructive measurement. Thus, 
the only remaining degree of freedom, within the chosen 
microcanonical ensemble, is the integer < K < 2S, that 
parameterizes Si = (K - S + N/2)/2 and S = (-K + 
S + N/2)/2 ineq.(||). The corresponding Gibbs distribu- 
tion over K with the temperature created by the cooling 
follows from the spectrum eq.(|^), which now takes a form 
E S mn(K) = (e + e/2)N + e(K - S) + of>/e) in the 
main e [i limit. It indicates that the ground state for 
given S, M, N corresponds to K = 0. The next excited 



K = 1 state within the microcanonical ensemble costs 
energy e. Thus, once temperature T < e, the micro- 
canonical ensemble for given S, M, N essentially consists 
of only one state — the ground state K = 0, given as 



\SMN) = ZsM N (al ] ) S+M (a^Y " (PJ) N/2 5 |0), (6) 

R = aoffflij, — ait<Joj.! 

where S = (S + N/2)/2 and Si = (-S + N/2)/2; 
Zsmn is the normalization coefficient, and |0) is the 
vacuum. These states realize the irreducible represen- 
tation of the isotopic SU (2) group. For S = 0, M — 
( TV even), the state (§) is the Sf7(2)-singlet. 

The state (|^) exhibits quantum fluctuations of individ- 
ual values of N aa , preserving, however, their three linear 
combinations: A 0T + N QL = S + N/2, A iT + A u = 
-S + N/2, N n +N ^-N u - N ol = 2M. Thus, there is 
only one degree of freedom, which we choose to be Nq-[. 
Most transparently, these fluctuations are described in 
the basis of the Fock states |iVoT , A^j., A^f, N^) , which 
arc the fragmented condensates |l. Expanding ^) over 
the Fock basis, we get the probability to find a given 
value of iVof : 



S-M 



Pi 



SAIN 



(JV 0T ) 



Z 2 SMN N ^(N/2 + S-N ^\ 
(A 0T - S - M)\{N/2 + M- A 0T )! 



(7) 



To estimate the strength of the quantum fluctuations, 
consider SN 2 p where for any Q: Q = (SMN\Q\SMN) 

and SQ 2 = (SMN\{Q~Q) 2 \SMN). At S = 0, we have 



SMN 



(iVof) = 2/ A over the whole range of < A( 



Of 



< 



N/2: and 5N 2 ^ = N 2 /12. For 1 < S < TV, S - \M\ > 
1, the Gaussian approximation of (^) gives 



A 0T 



(S + M)(S + N/2)/2S , 

(S 2 - M 2 )((N/2) 2 - S 2 )/8S 3 . 



(8) 
(9) 



N, we obtain SN^ 



For typical values S ^ 

A 2 /(8S) ~ A 3 / 2 . In a direct analogy with the s = 1 
case fa] , these anomalous fluctuations are consequence of 
the conservation of the total (iso)spin M. 

The state (|^) is a quantum superposition - a 
Schrodinger Cat — of the phase separated "classical" 
BEC's, each BEC being characterized by well-defined 
spatial density distribution of the spin-up and spin-down 
components. The phase separation is due to the involve- 
ment of two single-particle modes with different occupa- 
tions of the spin- up/down components. 

Measurement of JVof produces a collapse of the Cat. 
The resulting states essentially depend on how strong 
the interaction between the measuring device and the 
system is. The most "accurate" measurement of iVof, 
with the absolute error less than one particle, would col- 
lapse the Cat into a Fock state |JVbf , Aox ,Nif,Nti), with 
the probability (pi) to obtain given Aqt . Such a measure- 
ment, however, implies that the interaction of the system 
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with the apparatus is much stronger than that due to the 
Hamiltonian (^). In this sense, such a measurement is de- 
structive not only to the supposed-to-be- fragile Quantum 
Cat, but also to its classical counterparts. 

More interesting is a less destructive measurement, 
which collapses Quantum Cat into some "classical" (in 
a particular sense) states. These classical states are the 
phase-coherent BEC states, which can be interpreted in 
terms of the classical vectors So, Si interacting anti- 
ferromagnetically by the term 2/oiSo-Si = 7oi[S 2 — Sf — 
S 2 ], Iqi > 0. Comparing the classical energy with the 
term in the square brackets in the full (quantum) expres- 
sion we see that the creation of such classical states 
costs only the energy SEc ~ Iqi[So + S%] ~ IoiN ~ ^. 
Here we took into account that So ~ Si ^ A" for 
S< A. 

The above-discussed classical BEC states are readily 
obtained within the Gross-Pitaevskii approach. As an 
instructive example, let us construct the classical ground- 
states. To this end we replace the bosonic fields opera- 
tors by the c- fields. Accordingly, in eqs. (||,fl|), we replace 
the a-operators by the c-number amplitudes, and find 
the minimum of energy at given N, S, M. This yields 
«o T = VS + N/2 cos(/3/2) exp(z(a + 7 )/2), a ol = 
-y /S + N/2 s in(/3/2) exp(i(a - 7 )/2), a 1T 
iy/-S + N/2 sin(/3/2) exp(i(ai + 7)/2), aii = 
iy/-S + N/2 cos(/3/2) exp(i(ai — r ))/2). Here cos/3 = 
M/S, (ai — a )/2 = — e(l + o(fx/e))t; a is the global 
phase, and 7 determines the x, y-projections of the to- 
tal spin ( S x — — S sin [3 cos 7) . Comparison of these with 
eq. (|^) shows that the classical occupation numbers coin- 
cide with the expectations for the quantum ones. 

The classical fields are non-stationary, which is a conse- 
quence of the energy difference between the states ipo and 
(fi. Hence, wherever </?o( x ) and </?i(x) have substantial 
overlapping, the spatial densities of the two components 
experience strong fluctuations, the total density remain- 
ing constant. These fluctuations are of purely classical 
nature and have nothing to do with the above-discussed 
quantum fluctuations of occupation numbers! It is rea- 
sonable to reveal these classical fluctuations directly in 
the quantum wavefunction (^). To this end we calculate 
fluctuations of the difference of the component densities, 
r/(x) = *|(x)* T (x)-^|(x)*j.(x) [Note that rj^O for 
\M\ -C S, that is the state (g) does not exhibit any phase 
separation on average]. For the fluctuations of r\ we get 

W = 4[<Po + <Pi- ^WiW^ + 2¥>oVi (N + 4)- 

-S(S + l) + 4SW(M-So7)]. (10) 

The two terms in ([l0]) correspond to the quantum (~ 
SN^) and classical (~ N 2 for S <C N) fluctuations of 
77, respectively. Though the second term is dominant in a 
general case, it vanishes in the trap center, where ipi = 0, 
and the fluctuations are of purely quantum nature. 



To propose a particular set of measurements that could 
reveal the statistics (|?|) , we note that the quantities Nq — 
Ea N o, = N/2 + S, JVi = Y.o N i° = N/2 - S, A7 = 
N 0] + iVi T = N/2 + M, N x = N ol + N n = N/2 - M do 
not fluctuate. Thus, the values S, M, N can be ex- 
tracted, in effect, without destroying the wave function. 
Once the set of values S, M, N is fixed, one can mea- 
sure AT |, by, e.g., observing the density of the spin-up 
component in the center of the trap by the spatially 
differentiated selective imaging jl(|. In the trap cen- 
ter tpi — 0, therefore the component densities are un- 
ambiguously related to the occupations of the state <po 
|]ll| . Actual measurements always have some uncertain- 
ties AM, AS, AiV in determining AT, S, N, respec- 
tively. These may wash out the quantum fluctuations. 
First, we note that the mean value Nq^ in eq.(^) is 
exactly the value given within the classical geometri- 
cal picture. Thus, ((ATVot) 2 ), where (...) denotes 
the statistical mean, must be significantly smaller than 
the quantum fluctuation SN^p eq.@. Practically, for 
({AN) 2 ) < N, (S) ~ y/N, we find that, as long as 
((AM) 2 ) -c (S)/2, the quantum fluctuations can be well 
distinguished. 

We consider an intriguing alternative to the direct cool- 
ing of the two-component system in order to obtain the 
state ([|). Namely, the merging of two spatially separated 
condensed species from two identical wells with equal 
numbers of oppositely "polarized" particles in each of 
them. At T <C £ (which means, for our purposes, prac- 
tically zero), and in the absence of tunneling between 
the wells, the system is in its lowest energy state. We 
note further that our state is (N/2 + l)-fold degener- 
ate within the given sector of N and M = 0, the other 
N/2 groundstates being obtained by exchanging spin- up- 
spin-down pairs between the wells. The key point is that 
the initial state is a superposition of the eigenstates of 
the total pseudo-spin operator with the typical values 
< S < VN, and the degeneracy is exhausted by the 
isospin index (which, at M — 0, has exactly N/2 + 1 val- 
ues, from to N/2). Thus, each of the S-eigenstates in 
our superposition is not degenerate. Accordingly, when 
tunneling is slowly turned on and, finally, the inter-well 
barrier is adiabatically removed p"2[ | , the system remains 
in the superposition of the ground S-cigenstates, which 
at e/fi ^> 1 take on the form (0) ( M — 0). Projection 
on the particular S can be done by, say, measuring total 
density in the trap center jO]. 

Other (more exotic) options for the direct Cat-state 
engineering rely on merging more than two BEC clouds in 
the presence of the controlled SU (2) symmetry violating 
potential. We will analyze them elsewhere. 

Let us discuss the effect of the symmetry breaking static 
perturbation, which was shown to destroy the quantum 
fluctuating state in the case s = 1 in the thermo- 
dynamical limit 0. A similar situation occurs in the 
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case of the pseudo-spin s = 1/2 bosons. A role of 
the perturbation is played by the gradient B' of the 
pseudo-magnetic field B' = (0,0, B'), consisting of the 
difference U' = f/j — of the trapping potentials 
for the corresponding components and of the interaction 
term H = (g'/2) J dx [*{*|* T * T - *j*t where 
g' = ina'/m and \a'\ < a in 87 Rb §. 
model, we find the corresponding term 



Within the 



Hb — —B'(Sq z — Si z ), 

B' = Jd* (C/'(x) + 9 -p(*))(<p\ - ifl), (11) 

to be added to H' , eq.(|4|), where we take into account 
N « N, « AT/2, and p(x) = A(<^ 2 + ^)/ 2 . 

The term i/g does not commute with the total spin, 
which implies that it does not favor quantum fluctua- 
tions. It decreases energy of the classical (not fluctu- 
ating) state \BEC) = (aJt ^)^ 72 ! ) b y SE B bc ~ 
-(BEC\H B \BEC) w B'JV, ( B' > 0). At the same 
time, the classical states loose energy with respect to 
the lowest eigenstate of H' , eq.(Q), by the chemical 
potential ji. Accordingly, in order to insure the sta- 
bility of the Schrodinger Cat state the condition 
/i > |-B'|iV should hold. This condition cannot be obvi- 
ously satisfied in the thermodynamical limit, and should 
be considered in mesoscopic situations only. We note, 
however, that it can be well controlled within the cur- 
rent experimental capabilities. Taking the effective po- 
tential breaking the symmetry as the oscillator poten- 
mv 2 x 2 /2, with some frequency 
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tial SU = mv"x" I'l, with some frequency v 
find B' = v 2 /(2e) in eq.(ll), where we have employed 
the lowest eigenstates <fio : i, of the ID oscillator with fre- 
quency e. Accordingly, v should obey v/e< wJi/JeN). 
For N < 10 4 , (i < u>o, this yields v/e < 10~ 2 , which is 
well within the accuracy of the experiment |13|] . In fact, 
the control of the value B' together with the exter- 
nal rf-pulses [Q,[l(| provides a tool for manipulating the 
parameters of the Cat state. 

Another effect of the symmetry breaking, existing even 
for B' = 0, is opening the channel for irreversible tran- 
sitions 1 — > 0. While bringing the system to the abso- 
lute ground state ( S = N/2), it is, however, of a sec- 
ond order in g', and needs the normal component to 
be present to proceed at noticeable rate. In 87 Rb 0, 
g/g' w 30. Thus, a time scale for a such relaxation is a 
factor of 10 3 longer than relaxation time ~ 1/g 2 of any 
lowest energy excitation at considered low temperatures 
and densities, which makes it irrelevant for all practical 
purposes, if compared with the time of the Cat-state for- 
mation ~ 1/fi ~ 1/g . 

Let us consider stability of the above analysis with re- 
spect to external symmetry breaking temporal noise. For 
example, the trapping potential may be slightly different 
for the components and can fluctuate in time. These fluc- 
tuations may produce two qualitatively different effects: 



1) heating of the cloud at some rate r^ 1 , and 2) dis- 
ruption of the coherence at a rate Ty r , which, generally 
speaking, is different from t// 1 . 

A most detrimental for the Cat situation may occur 
when the noise, while producing no significant heating, 
disrupts the anomalous quantum correlations A(t) = 

As dis- 



SNqp as given by eq 
cussed in ref. 



at times i~; rr <C 7V 



|l4f , this occurs when the noise correla- 
tion length is larger than a typical thermal length of the 
cloud. Accordingly, if the noise is produced by small cor- 
relation length factors, such as, e.g., collisions with the 
background gas, the decoherence is not faster than the 
heating. 

In order to investigate the case under consideration, 
we note 
that A(t) can naturally be represented in terms of the 
one-particle pW(x,x',i) = (*t( X) i)$( x ', t)) (OPDM) 
and two-particle p^^xi, x 2 , x.' 2 , x' l7 1) = 

(* t (xi,t)^t( X2j t)^(x^,t)*(x' 1 ,t)) (TPDM) density 
matrices, respectively, where the averaging is performed 
over the initial state, which is taken to be the Cat-state 
(^), and over the noise. Accordingly, the equations for 
them should be analyzed and the rates t~*, t// 1 eval- 
uated and, then, compared with each other. In the case 
of the white noise, this can be done, practically, exactly 
with the help of the Furutzu-Novikov theorem as dis- 
cussed in ref. fll4| . 

The heating rate is defined as the inverse time Th re- 
quired to increase the total kinetic energy per particle 
K(t) = -(l/m)V2 =x , (0 ( 1 )(x,x / ,i)/A ~ t by e. The 
decoherence rate r" 1 is defined as the inverse of a typ- 
ical shortest time of the decay of the elements of the 
OPDM and the TPDM. All rates can be evaluated with 
respect to the OPDM and the TPDM solutions in the 
limit t -> 0. 

It is important to note that, while the heating is most 
effectively produced by a short range noise, the coherence 
is destroyed by the long range noise due to the long range 
temporal fluctuations of the trapping potential U'(x, t). 
These fluctuations can be produced by thermal or any 
other noise of electric currents I'(t) in the magnetic 
coils. For the purpose of estimating, we assume that 

U' = U(x)l'(t)/I, (I'(t)l'(t')) = C 2 6(t - t'), (12) 

where U (x) is some typical trapping potential, which 
does not fluctuate; I stands for the static current; we 
assumed the white noise structure of the current fluctu- 
ations, which is characterized by the constant C 2 . In 
the case of large temperatures of the coils Tq, this con- 
stant is given by the fluctuation-dissipation theorem as 
C 2 = 2Tc/R, where R stands for the resistance. Pro- 
ceeding as described above and in ref. juj and taking 
into account eq.jl 7 ^), the heating and the decoherence 
rates can be found as 



(V[/(x)) 2 C7(meI 2 ), 



U 2 (x)C 2 /I 2 , (13) 
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ature), e/h ss 10 J s , we find r h 
if the coil power RI 2 
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-15 r 



(e/h) 2 r c (I 2 )/I 
" 6 s, e/7i « 10 3 s" 



respectively, where the line stands for a typical average 
over the cloud volume, and the numerical coefficients of 
order 1-10 (depending on the order of the correlator) are 
omitted. 

We estimate U 2 {x) « e 2 and take into account that 
typical kinetic energy is of the order of the potential 
energy, that is, (VC/(x)) 2 /to ps e 3 . For the thermal 
noise, where C 2 = 2T C /R, T c = 300K (room temper- 

' pa e 2 k B T/{h 2 RI 2 ) ~ 
1W (we restored the 
standard units). Accordingly, the thermal noise of the 
coil currents can be ignored for all practical purposes. 

In a general case of a non-thermal noise characterized 
by some relaxation time tq and the current fluctuation 
magnitude (i' 2 ) , so that (/'(*}/'(*')) = (i' 2 ) exp(-|t - 
t'\/ T c), the white noise result juj and the estimate ( |l3| ) 
can still be employed as long as tc -C Ti rr . We, then, 
take C 2 ps t c (I 2 ) , and obtain ps r^/%\2^/ T '2\ 
This implies that, for , e.g., tc = 10 
the relative fluctuations (I 2 ) /I 2 can be as large as 10% 
and yet n rr > 100s, which should be compared with the 
time of the Cat-state formation |Q. As discussed above, 
this time can be as short as ~ 1/fi ~ 10~ 2 — 10~ 3 s. 

The decoherence due to losses can similarly be ana- 
lyzed within the framework of equations for the OPDM 
and TPDM. Regardless of the nature of the interaction 
vertex responsible for the losses, these equations can be 
obtained as an expansion with respect to the vertex, and, 
then, the rate t~1 compared with the corresponding 
contributions to the rate of losses Tj", , defined as the 
time derivative of J dxpW (x, x, t) per one particle. Per- 
forming this general analysis, we find that t~* and 
T ioss are °f the same order. Thus, once the losses are 
kept low during the time of the experiment (the scale 
of which is determined by the adiabaticity requirement 
IP^I ), the decoherence can be ignored. 

In summary, we have discussed the role of the isotopic 
SU(2) symmetry in the formation of the ground states 
in a confined mesoscopic two-component Bose gas. The 
conservation of the total isotopic spin results in conden- 
sation of a thermal cloud, consisting of the two uncorre- 
cted components, on two lowest single particle states in 
the trap. The resulting ground state is essentially non- 
classical Schrodinger Cat state characterized by anoma- 
lously large non-uniform fluctuations of the density of 
each component about its classical expectation. The 
symmetry breaking term may destabilize this state into 
the phase separation, with both components becoming 
the classical BEC. This process, however, can be well 
controlled experimentally The realization of the fluctu- 
ating many body state provides unique opportunity for 
observing the effect of measurement induced collapse of 
the Schrodinger Cat to classical states. The conditions 
required for distinguishing the quantum fluctuation from 
the statistical noise are delineated. 
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